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The local distribution of exciton levels in disordered cyanine-dye-based molecular nano-aggregates 
has been elucidated using fluorescence line narrowing spectroscopy. The observation of a Wigner- 
Dyson-type level spacing distribution provides direct evidence of the existence of level repulsion of 
strongly overlapping states in the molecular wires, which is important for the understanding of the 
level statistics, and therefore the 'functional properties, of a large variety of nano-confined systems. 
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One of the current dreams in the field of molecular op- 
tics is the full understanding of nature's way to harvest 
and use photonic energy which ultimately could enable 
the development and design of highly efficient functional 
optical devices using molecular arrangements as build- 
ing blocks. One of the crucial elements of such photonic 
assemblies are the 'wires' which transport the energy be- 
tween the different functional units of the devices. Natu- 
ral systems often utilize structures of coupled aggregated 
pigments to transport energy in the form of excitonic ex- 
citationSfi"— Such structures can also be mimicked in syn- 
thetic systems, greatly assisting studies aiming to under- 
stand their fundamental properties. An important class 
of synthetic species, on which we focus here, is found in 
the so-called one dimensional (ID) J-aggregates based 
on, for instance, pseudoisocyanine, porphyrin, and benz- 
imidazole carbocyanine dyesi^"— 

Synthetic, as well as natural systems, usually exhibit 
a substantial degree of disorder, arising from the envi- 
ronment and from vibrations and disorder within the 
systems themselves. In general the presence of disor- 
der in gapped systems leads to the formation of highly 
localized states inside the optical, electronic or mag- 
netic energy gap of the unperturbed system; i.e. to 
a tail of the density of states inside the gap generally 
referred to as the Lifshits tail4 There are many sys- 
tems which optical properties are governed by exciton- 
like excitations highly susceptible to disorder leading to 
localization and level repulsion phenomena. These in- 
clude conjugated oligomer aggregates^ and polymers,— 
molecular J-aggregates,—"— semiconductor quantum wells 
and quantum dots,— gold nanoparticles,— semiconduc- 
tor quantum wiresj^ as well as photosynthetic light har- 
vesting complexeaii^ and proteins^ (see Ref . for a recent 
overview). In all these systems, excitons are confined at 
least in one dimension at a nanometer scale. 

The physical and transport properties of most of the 
above mentioned systems are predominantly determined 
by the states residing in the vicinity of the energy gap, 
i.e. the gap excitation itself and the Lifshits tail be- 
low it, even at finite temperatures The localization of 
the exciton states within the Lifshits tail gives rise to a 



local (hidden) statistics of the levels, which deviate sub- 
stantially from the overall statistics.— Therefore, un- 
derstanding the physical properties of these systems re- 
quires the use of statistical approaches; the energy level 
distributions become an important part of the theory and 
interpretation of the experimental data. 

For non-interacting systems it is well known that the 
presence of disorder leads to an energy spectrum with 
a Poissonian energy spacing distribution. In less triv- 
ial cases of interacting systems, the situation naturally 
becomes more complex leading to the concept of level re- 
pulsion, i.e. a vanishing probability to find two quantum 
states with the same energy. The level statistics of such 
a system is known as Wigner-Dyson statistics (see the 
excellent textbook by Metha Ref. [13 for an overview) . 

Level repulsion phenomena in nano-confined materi- 
als has recently drawn considerable attention, in par- 
ticular, concerning localized Wannier excitons in disor- 
dered quantum wel b^^d^ and wires;^ and in disordered 
graphene quantum dotS(^ as well as concerning vibronic 
states in polyatomic molecules;^ Time-resolved resonant 
Rayleigh scatteringi^ and near-field spectroscop y^^'^" 
have been used to study them. In Ref. |23| . an alternative 
method has been proposed to analyze the level statistics 
- low-temperature time-resolved selectively excited ex- 
citon fluorescence spectroscopy, widely known as fluores- 
cence line narrowing (FLN) spectoscopy. Under a narrow 
(compared to the the J-band width) excitation within 
the J-band, the fluorescence spectrum consist of a sharp 
intensive peak at the excitation energy and, growing in 
time, a red-shifted feature, resulting from the exciton 
band relaxation. It is this feature that contains informa- 
tion about the level statistics of the spatially overlapped 
states. 

Here we apply a variant of this method to reveal the 
repulsive statistics of levels residing in the Lifshits tail^ of 
J-aggregates of pseudoisocyanine (PIC) with a chloride 
counter- ion (PIC-Cl). In contrast to the earlier proposal 
to use time dependent FLN^^ we here show that also 
steady state FLN, which is much more simple in real- 
ization, can be utilized to extract the desired informa- 
tion. We extract the conditional probability distribution 
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of the repelled states from the experimental data and 
show that this distribution is Wigner-Dyson-like indicat- 
ing zero probability for zero energy spacing, which is a 
fingerprint of the level repulsion. 

Figure [T] shows the absorption A and fluorescence F 
spectra of J-aggregates of PIC-Cl, the latter measured 
after excitation using 400 nm light. Both spectra exhibit 
an intense peak arising from the dominant exciton tran- 
sitions (the J-band) and a much less intense and broad 
shoulder located on the red side of the J-band. We relate 
this red feature to aggregates in the vicinity of the sub- 
strate as the relative intensity of this shoulder decreases 
upon increasing thickness of the aggregated film. For 
completeness, we note that the overall absorption spec- 
trum of our samples (see inset in Fig. [T]) is found to be 
in good agreement with earlier resultsi^ 
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FIG. 1: Low temperature steady-state absorption (solid line) 
and fluorescence (dotted line) spectra of J-aggregate of PIC- 
Cl in the neighborhood of the J-band. The fluorescence spec- 
trum was measured after off-resonance excitation far in the 
blue tail at temperature T — 4K. The inset show the absorp- 
tion spectrum in a wider spectral range. 

Before turning to the main experimental results, we 
briefly sketch some of the theoretical background relevant 
for the present study; more details on the model and the 
energy level structure of the disordered J-aggregates has 
been discussed in Ref. [l^. A typical realization of the 
calculated low-energy level structure and wavefunctions 
for a one dimensional aggregate of 300 chromophors is 
depicted in Fig. [2] In this calculation, we used a gaus- 
sian disorder distribution of the chromophore energies 
with standard deviation a = O.IJ (disorder degree from 
now on), J being the transfer interaction between chro- 
mophores (for more details see Ref.[ii). 

Without disorder wavefunctions are fully delocalized, 
and the lowest state is located at -2.404J. The pres- 
ence of disorder leads to localization of the wavefunc- 
tions within so called segments, and to the appearance 
of highly localized states within the band gap, i.e. in 
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FIG. 2: A typical realization of the exciton wave functions ipun 
{v — 1 . . . 14) in the neighborhood of the bare exciton band 
edge Eb/.J = -2.404. The Lifshits tafl of the DOS {E < Et) 
is shaded. The origin of the energy is chosen at Eq = 0, the 
baseline of each wavefunction represents its energy in units 
of J. Wave functions are in arbitrary units. N* denotes the 
typical localization size of the tail states. Filled red curves 
are s-like states which overlap weakly. Some of them appear 
even slightly above the bare band edge. Filled blue curves 
are p-like states which overlap well with their s-like partner 
state lying below. Higher grey-shaded states are band states. 
They are delocalized to a larger extent as compared to the 
tail states. 



the Lifshits tail (grey shaded area in Fig. [2]). These are 
the states of our primary interest since they determine 
the optical properties and transport in the J-aggregates. 
They originate from localization in well-like fluctuations 
of the site potential on the molecules. The optically dom- 
inant states resemble s-like wavefunctions which have 
no nodes within their localization segments. The s-like 
states lying deep in the Lifshits tail usually appear as sin- 
glets and are localized by the so called optimal fluctua- 
tions of the site energy)^ Close to the band edge, however, 
the s-like states often have partners localized within the 
same localization segment. The latter look like p-states, 
having one node within their localization segment. Man- 
ifolds like these form the local (hidden) structure of the 
tail of the density of states. Since these states are local- 
ized on the same segment, one may expect level repulsion 
to occur for them, as is indeed observed. In contrast, 
states from distant (non-overlapping) manifolds can be 
arbitrarily close in energy. 

Optical experiments probe the states with a finite tran- 
sition dipole moment. For the s-like states, the transi- 
tion dipole moment scales proportionally to V N* , where 
N* is the typical localization length of the states. This 
enhancement of the dipole moment is known as superra- 
diant enhancementi^ Typically, the p-like states have a 
transition dipole moment which is several times smaller 
than the s-like statesi^ Nevertheless, since the p-states 
are not perfectly antisymmetric, they do have a finite 
transition dipole moment and these states can be opti- 
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cally excited too. Therefore optical experiments can be 
used to probe the level statistics by studying the relax- 
ation between p and s like levels.— 
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FIG. 3: Low-temperature (T = 4 K) fluorescence spectra 
of J-aggregates of PIC-Cl measured while selectively exciting 
within the J-band. The vertical lines show the positions of the 
red feature maxima. The excitation wavelengths are indicated 
along the right vertical axis. 

To experimentally study the level statistics in the 
neighborhood of the exciton band edge, we performed 
steady state resonance fluorescence measurements using 
a narrow excitation line to excite states within the J- 
band. Figure [3] shows a number of such spectra recorded 
at low temperature using different excitation energies. 
The spectra show a strong peak at the excitation wave- 
length, together with a broad red-shifted emission band 
separated from the main peak by a pronounced dip. The 
red shifted emission originates from relaxation of the ini- 
tially excited exciton states into states of the Lifshits 
tail of the DOS. The position of the maximum of this 
band remains almost unchanged while J-aggregates are 
excited on the blue side of the J-band. For the red- 
side excitation, the peak position moves to the red, the 
line shape changes considerably, and the dip washes out. 
This is a consequence of the changes in the relaxation 
pathways since red-side excitation predominantly excites 
s-like states. 

The dip close to the excitation energy in blue side ex- 
cited spectra distribution shows that energy relaxation 
into states close to the excitation energy is substantially 
suppressed, hinting to the occurrence of level repulsion. 
The sheer existence of the dip, however, is not enough to 
conclude on the level statistics. The problem is that the 
relaxation process from the initially excited states to the 
lower lying levels is phonon assisted and, hence, the line 
shape is determined by the product of the level spacing 
distribution function and the phonon spectral density. 
Since the phonon spectral density vanishes for zero en- 
ergy, one expects the spectral intensity to vanish at the 
excitation energy, even without level repulsion. More- 
over, one also should bear in mind that together with the 



fluorescence of the relaxed excitons two more processes 
contribute to the red feature and affect its line shape: the 
phonon side-band fluorescence?^ and surface-mediated 
fluorescence in our thin samples. All three contributions 
to the red-shifted feature are spectrally superposed and 
must be separated in order to extract the signal we are 
interested in. We note, however, that if the observed red- 
shifted feature would solely originate from the phonon 
side-band, its line shape and position would be virtually 
independent of the excitation wavelength, clearly in con- 
tradiction to the experimental observations. 

In the analysis of the observed spectra we limit our- 
selves to those spectra measured using excitation on the 
blue side of the J-band, since it is here that one ex- 
pects the p-like states to contribute most strongly. In 
order to discriminate the true relaxation mediated flu- 
orescence (RMF) from the surface mediated and the 
phonon sideband fluorescence we use a simple subtraction 
method. For this, we consider the differential spectrum 
between two experimental spectra with close excitation 
wavelengths A2 > Ai defined by 

AF(Ai,A2,A) = F(A2,A)-/3F(Ai,A-A2 + Ai) , (1) 

where the second term on the right hand side is the 
F(Ai, A) spectrum shifted in wavelength to match its ex- 
citation peak position with that of the J"(A2, A) spectrum. 
In addition, this term is rescaled by a factor of (3 in or- 
der to cancel the red tail in the spectra; any feature that 
is not wavelength dependent is suppressed in the differ- 
ence spectrum [U and the resulting difference spectrum 
represents just the RMF differential signal AR{Xi, A2, A): 



Ai?(Ai,A2,A) = i?(A2,A)-/3i?(Ai,A) , 



(2) 



At the next step, we calculated the quantum efficiency 
of the red-shifted feature and found that it did not ex- 
ceed 0.3 for most blue excitation. For the spectra we will 
use in the fitting procedure, the efficiency even smaller, 
around 0.1. This means that excitons make only one 
step of relaxation, moreover, the major contribution to 
this process comes from the intra-segment hops (see the 
discussion in Section 2). Then the theoretical RMF line 
shape i?(Ae, A) ^ S{X — Ae)-Psp(Ae, A — Ag), and we can 
relate two RMF spectra taken for different (close) exci- 
tation wave length A2 > Ai as: 



i?(A2,A) 



5(A- A2) 



i?(Ai,A) 



(3) 



^(A-Ai) 

where we assumed that the energy spacing distribution 
function varies much slower than the phonon spectral 
density, the assumption which, as will be seen, is consis- 
tent with the final results. Substituting Eq. ([3]) into ©, 
we arrive at a relationship between the differential and 
ordinary RMF spectra: 



where 



Ai?(Ai,A2,A)«5(A)i?(A2,A) 
S{X-Xi) 



?(A) = 1-/? 



^(A- A2) 



(4) 
(5) 
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i.e., the lineshape of the RMF spectrum can be extracted 
from the Hneshape of the differential RMF spectrum by 
dividing it by the known correction function, provided 
the factor fi is adjusted to cancel long red tale. 

It is important that the correction function is almost 
constant for wave lengths that are far from the excitation 
wave lengths (A — Ai,2 ^ IA2 — Ai|). Hence, it would not 
change the shape of the distant features in the long red 
tail of the experimental fluorescence spectra. Such fea- 
tures will therefore be canceled in the difference spectrum 
which will contain only the contribution of the RMF. Fi- 
nally, applying the above reasoning and formulae to the 
experimental FLN spectra, we can recover the MRF line- 
shape i?(Ae, A) according to the formula (jH). 



for details Ref. |23 
presented in Fig. [TJ 



which reproduces the Stokes shift 
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FIG. 4: The conditional probability of the nearest level spac- 
ing distribution Psp (diamonds) obtained by dividing the ex- 
perimental curve (circles) by the Debye-like spectral density 
S(A) oc A^^ together with the calculated level spacing dis- 
tribution Psp (dashed curve). The experimental RMF spec- 
trum (circles) show the data after applying the subtraction 
method (described in the text) to eliminate the contribution 
of the long red tail resulting from non-RMF transitions. Also 
shown are the theoretical RMF spectrum calculated for the 
Debye-like spectral density ^(A) oc (solid curve) and the 
absorption spectrum (dotted curve). 

The described procedure has been applied to two 
spectra recorded using excitation in the blue part of 
the absorption spectrum (Ai = 568.5 nm and A2 = 
569 nm). The result, assuming a Debye-like spectral den- 
sity S{\) (X A~^ is shown in Fig. H) The extracted RMF 
spectrum, obtained using (3 = 1.06, is shown in the figure 
by the open circles. Clearly, the non-RMF contributions, 
leading to the long red tail the fluorescence spectrum, is, 
as expected, nearly fully eliminated. Superimposed on 
this experimental spectrum is the result of a simulation 
of the RMF fluorescence spectrum, using a Gaussian dis- 
order with standard deviation a = 0.2, again assuming 
a Debye model for the phonon spectral density, and an 
exciton-phonon scattering strength of Wq = 22.4 J (see 




FIG. 5: Open circles represent the experimental red-shifted 
feature obtained after applying the subtraction approach (de- 
scribed in the text) to eliminate a contribution of the long red 
tail resulting from non-RMF transitions. The solid curve is 
the theoretical RMF spectrum calculated for the Debye spec- 
tral density 5(A) oc A^'^, while the dashed-dotted and dashed 
curves are the RMF spectra for S(A) oc A and S{\) oc A~^, 
respectively. The dotted curve denotes the absorption spec- 
trum. 

It is to be noticed that calculations performed with 
non-Dcbye models for the phonon spectral density, do 
not lead to a satisfactory agreement with the experimen- 
tal data. This is clearly demonstrated in Fig. [51 The 
apparent validity of the usage of the Debye model cor- 
roborates the results of Ref. [13, where this model has 
been successfully used to explain the temperature de- 
pendence of the J-band width and the radiative lifetime 
of J-agregates of the dye pseudoisocynine wiith different 
counter ions. 

Figure m shows the principal result of the present work: 
the conditional distribution of the nearest-level spacing, 
^sp(Ae,A — Ae) (diamonds), obtained after dividing the 
extracted RMF spectrum presented in Fig. 2] (open cir- 
cles) by the phonon spectral density S{\) ~ A""^. We see 
that Psp(Ae,AA) tends to zero upon AA — > 0. This is 
a clear signature of the repulsive statistics of the near- 
est level spacing, i.e., it is of a Wigner-Dyson-type. In 
spite of the observed distribution is strikingly close to 
a Wigncr-Dyson one, it remains difficult to determine 
whether it is the truUy Wigner-Dyson distribution, re- 
quiring a linear decrease to zero probability for zero spac- 
ing, which, though consistent with the data, is not fully 
proven by the current experiments and analysis. 

To conclude, we experimentally studied the statis- 
tics of the low energy spectrum of disordered molecu- 
lar nano aggregates of pseudoisocyanine with the chlo- 
ride counter ion in the neighborhood of the exciton band 
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edge. The fluorescence line narrowing technique, allow- 
ing to probe the local energy level distribution;^ has 
been exploited for this goal. We found a clear signature 
of a Wigner-Dyson-like distribution for the nearest level 
spacing, originating from the exciton states localized on 
the same segment of the aggregate and thus undergo- 
ing the quantum mechanical level repulsion. This is the 
first direct experimental prove of the existence of hidden 
structure of the exciton low energy spectrum, the region 
which dominates the aggregate optical response and low- 
temperature transport. 

Finally, we note that our finding has a wider applica- 
bility than for a simple ID Frenkel exciton system consid- 
ered here. The reason is that the nature of the band edge 



states (mostly in the Lifshits tail) is shared by a large va- 
riety of systems, such as gold nanoparticlcsi^, quantum 
wells and quantum wires^Sr^. 
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